clear;clc;
model = load('T_S_Wr.dat');
T = model(:,1);
S = model(:,2);
Wr = model(:,3);
% plot(T,S)
plot(T,Wr);
Pz = load('Pz.dat');
XZ = load('X_Z.dat');
X = XZ(:,1);
Z = XZ(:,2);
wigb(Pz,X,Z,1);